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Abstract. 


An  acoustic  Doppler  current  profiler  (ADCP)  was  applied  to  measure  the 
three-dimensional  flow  field  in  a  tidal  channel  along  the  coast  of  Vancouver 
Island,  British  Columbia.  The  ADCP  was  rigidly  mounted  to  the  bottom 
in  the  center  of  the  channel,  and  rapidly  sampled  the  velocity  profiles  along 
its  inclined  beams.  The  data  were  directly  read  by  a  computer  on  shore. 
This  method  of  deployment  allows  for /an  explicit  estimate  of  the  uncertainty 
of  mean  velocity  due  to  turbulence,  and  obtains  the  weak  signal  of  vertical 
velocity.  The  mean  flow  vector  is  derived  from  the  measured  beam  velocities 
under  the  assumption  that  the  flow  is  statistically  homogeneous  in  the  hor¬ 
izontal  plane  over  the  distance  separating  the  inclined  beams  of  the  ADCP. 
It  is  argued  that  sufficient  averaging  is  required  for  flow  measurements  in  a 
highly  turbulent  environment. 

The  depth-time  variations  of  velocity  and  shear  in  this  natural  channel 
are  more  complex  than  is  found  in  idealized  one-dimensional  channel  flow, 
which  seldom  occurs  in  nature.  The  channel  contains  strong  secondary  cir¬ 
culation,  intense  up-  and  down-welling  intervals,  frequent  shear  reversals  and 
a  substantial  amount  of  transverse  velocity  shear. 
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1.  Introduction 


Flow  in  shallow,  coastal  waters  is  complicated  in  its  space-time  variations. 
Good  current  measurements  are  essential  for  understanding  coastal  dynam¬ 
ics,  developing  and  testing  numerical  models,  and  solving  environmental 
problems.  The  most  significant  flow  component  is  usually  the  tidal  current, 
which  is  turbulent,  and  is  influenced  by  the  complex  bottom  topography 
and  coastal  geometry  typically  found  in  natural  channels.  For  example,  ob¬ 
servations  have  shown  that  transverse  circulation  is  generated  by  the  flow 
curvature,  and  flow  separation  can  occur  near  headlands  (Geyer,  1993)  or  in 
the  lee  of  islands  (Wolanski  et  al.,  1984).  The  transverse  flow  and  flow  sepa¬ 
ration  may  induce  up/down- welling,  by  the  mechanisms  proposed  by  Garrett 
and  Loucks  (1976),  and  Wolanski  et  al.  (1984). 

Oceanic  current  are  conventionally  measured  by  deploying  various  types 
of  current  meters.  A  current  meter  measures  the  flow  at  a  fixed  position, 
hence  an  array  of  meters  are  needed  to  study  the  flow  structure.  The  advent 
of  acoustic  Doppler  current  profilers  (ADCP)  provides  an  efficient  way  to 
measure  velocity  profiles.  The  spatial  resolution  and  range  in  profiling  are 
usually  adequate  to  measure  the  flow  and  shear  throughout  a  large  portion  of 
the  water  column  in  coastal  seas.  However,  the  velocity  measured  by  a  stan¬ 
dard  ADCP  is  the  radial  speed  along  each  of  its  four  beams  which  are  inclined 
30°  to  the  centerline  of  the  instrument.  The  true  velocity  vector  is  derived 
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from  the  measured  velocities  along  beam  pairs.  The  derivation  assumes  that 
the  flow  is  homogeneous  in  horizontal  plane  over  the  distance  separating  the 
beams  (e.g.,  Lohrmann  et  al,  1990).  In  a  turbulent  environment,  this  as¬ 
sumption  may  only  hold  statistically.  Hence  averaging  is  required,  and  only 
the  time  mean  velocity  vector  can  be  obtained.  From  this  point  of  view,  an 
AD  CP  is  not  equivalent  to  an  array  of  current  meters  equally  spaced  over  the 
profiling  range,  as  the  latter  does  not  involve  the  homogeneity  requirement. 
Turbulence  also  introduces  a  statistical  uncertainty  in  the  mean  velocity  esti¬ 
mates,  hence  each  mean  ensemble  must  contain  sufficient  samples  to  reduce 
the  uncertainty  to  an  acceptable  level.  The  vertical  flow  is  usually  weak, 
requiring  that  the  ADCP  is  mounted  rigidly  so  that  this  component  is  not 
influenced  by  the  platform  motions. 

The  results  to  be  presented  are  based  on  data  collected  with  an  ADCP 
mounted  rigidly  to  the  bottom  of  Cordova  Channel,  British  Columbia.  This 
contribution  describes  the  “mean”  flow  in  the  channel,  i.e.,  those  with  time 
scales  long  compared  to  the  turbulence.  Section  2  describes  the  experiment 
setup.  Section  3  provides  the  calculation  algorithm.  In  section  4,  we  present 
a  detailed  data  analysis,  to  determine  the  uncertainties  in  the  mean  velocity 
estimates  from  various  sources,  and  to  justify  the  averaging  time  chosen  for 
the  calculation  of  the  mean.  The  observed  depth-time  variation  of  the  flow 
and  shear  are  presented  in  section  5.  The  deployment  methods  and  the 


3 


channel  flow  dynamics  are  discussed  in  section  6.  A  summary  is  provided  in 
section  7. 

2.  Experiment  setup 

The  measurements  were  taken  in  Cordova  Channel  in  September  19-30,  1994, 
during  a  multi-investigator  experiment.  A  broadband  600  KHz  ADCP  was 
located  in  the  center  of  the  channel  where  it  is  about  1  km  wide  and  30  m 
deep  (figure  1).  The  coast  line  at  the  western  side  of  the  channel  presents  a 
headland  (Cordova  Spit),  with  shallow  water  in  Saanichton  Bay  to  the  north 
(see  the  20-m  isobath).  At  the  eastern  boundary  of  the  channel,  James  Island 
provides  a  small  curvature,  and  both  the  northern  and  southern  tips  of  the 
Island  are  within  the  tidal  excursion  distance  from  the  site  of  ADCP.  The 
water  depth  is  fairly  uniform  along  the  axis  of  the  channel. 

The  transducers  on  the  ADCP  were  held  1.59  m  above  the  bottom  by 
a  metal  quadri-pod  sitting  on  the  seabed.  The  tilt  sensors  of  the  ADCP 
measured  heading,  pitch  and  roll  angles  with  0.01°  resolution  and  the  readings 
were  steady  throughout  the  experiment.  Hence  the  instrument  was  rigidly 
mounted.  Because  the  magnetic  field  was  distorted  by  the  metal  stand,  the 
compass  readings  were  not  correct,  the  true  heading  angle  was  determined 
to  be  150.3°  through  a  post-experimental  calibration.  A  shore  cable  800 
m  long  connected  the  ADCP  to  the  power  supply  and  a  computer  in  the 
shore  station,  and  was  used  to  send  instructions  and  to  communicate  data. 
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A  total  length  of  4.5  days  of  data  were  collected  using  the  routine  working 
mode  (mode  4)  of  the  ADCP  which  provides  the  maximum  profiling  range. 
The  profiling  range  was  limited  to  27.6  meters  above  the  bottom  by  the  side 
lobes  from  the  transducer  reflected  from  the  surface.  The  profiling  range  was 
broken  up  into  uniform  segments  (called  velocity  bins)  with  1  m  spacing. 
The  position  of  the  center  of  the  first  velocity  bin  above  the  bottom  was 
restricted  by  the  height  of  the  quadri-pod  and  a  1.5  m  blanking  distance 
after  transmission.  The  centers  of  useful  bins  were  3.6-27.6  m  above  the 
bottom. 

The  velocity  was  recorded  in  beam  coordinates  with  each  beam  nominally 
inclined  by  30°  from  the  vertical.  The  ADCP  pinged  at  the  fastest  rate 
possible,  about  0.75  second  per  ping  or  at  a  frequency  of  1.3  Hz.  The  number 
of  pings  averaged  by  the  ADCP  into  ensembles  ranged  from  2  to  4.  The 
detailed  data  analysis  presented  in  section  4  are  based  on  the  4-ping  averaged 
ensembles  collected  between  day  23.9  to  25. 

Throughout  the  experiment,  winds  were  typically  less  than  3  m  s-1  and 
only  reached  5  m  s-1  occasionally,  hence  the  surface  stress  was  negligible. 
A  moored  CTD  observed  a  general  warming  trend  of  1°C  coupled  with  a 
decrease  in  salinity  of  1.5  psu  over  the  9  days  of  the  experiment,  and  also 
some  moderate  short-term  variability.  CTD  profiles  obtained  by  the  IOS 
(Institute  of  Ocean  Sciences,  Patricia  Bay,  B.C.)  group  show  that  the  strat- 
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ification  in  the  channel  changes  with  time,  and  the  water  below  mid-depth 
(about  15  meters)  is  fairly  homogeneous.  Currents  peak  at  about  1  m  s-1. 
During  the  experiment,  the  tide  changed  from  spring  to  neap  and  the  diurnal 
constituents  became  increasingly  dominant. 

3.  Calculation  algorithm 

The  velocity  along  a  beam  at  any  cell  is  the  vector  sum  of  the  horizontal  and 
vertical  components  at  the  position  of  that  cell.  The  beam  configuration  is 
illustrated  by  figure  2,  with  <^1,2,3  representing  the  heading,  pitch  and  roll 
angles  respectively.  Denote  the  velocity  along  the  i—  th  beam  as  bt  (i  =1  to 
4)  and  the  velocity  vector  at  the  position  of  6,-  as  (u;,  u,-,  wt).  (There  should 
be  one  more  index  denoting  the  order  of  bins,  or  depth  cells.  This  index 
is  dropped  because  we  are  considering  only  a  single  bin.)  At  first,  we  will 
not  consider  the  heading  angle  because  it  only  involves  rotation  in  the 
horizontal  plane.  Then  for  small  pitch  and  roll  angles  and  correct  to  the  first 
order  in  (p2  and  y?3,  one  obtains  (Lohrmann  et  al,  1990) 

hi  =  — ui(sin  6  +  (f3  cos  0)  —  tui(cos  0  —  <p  3  sin  0)  +  vi<p2  cos  0 
b2  =  U2(sin  0  —  (ps  cos  0)  —  w2( cos  0  +  (f  3  sin  0)  +  v2ip2  cos  0 
63  =  —  u3(sin#  —  (p2  cos  0)  —  w3(cos  0  +  <p2  sinf?)  —  U3<^3  cos  0  (1) 

64  =  u4(sin  0  +  <£2  cos  0)  —  w4(cos  0  —  <p2  sin  0)  —  u4ip3  cos  0 
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where  9  =  30°  is  the  beam  inclination  angle  with  respect  to  the  centerline  of 
the  ADCP. 


We  define  four  quantities  as 


„  _  i>2  —  b\  b\  +  b'2 

U  2  sin  9  2  cos  9 

a  b4  —  b3  b3  +  b4 

2  sin  6  2  cos  9 

A  b\  +  b2  +  b3  +  b4  b2  —  b\  b4  —  63 

w  = - : - 7 - <p3  0  +  V?2- 


4  cos  9 


e  = 


2  sin#  2  sin# 

(61  +  62)  —  (63  +  #4) 

4  cos# 


(2) 


where  (u,v,iu)  and  e  are  combinations  of  the  measured  beam  velocities. 
When  the  ADCP  operates  in  the  “earth  coordinates”  mode,  the  transfor¬ 
mation  in  (2)  and  a  compass  rotation  is  applied  to  the  echo  from  every  ping 
to  produce  the  “velocity  vector”  and  the  “error  velocity”.  We  operated  in 
“beam  coordinate”  mode  in  which  no  transformation  are  applied  and  the 
beam  velocities  are  returned  directly.  The  relationship  between  e) 

and  ( Ui ,  Vi,  W{),  correct  to  the  first  order  in  < p2  and  ip3,  can  be  derived  to  read 


u  — 


Ui  +  u2 

W2  —  Wi 

u2  -  Ui 

+  <^2 

2 

2  tan  9 

<l°3  sin  29 

v3  +  v4 

w4  —  w3 

,  «4  -  ^3 

-  ¥>3' 

2 

2  tan  9 

+  V,2sin2« 

w  = 


V2  ~  Vl 

2  tan  9 
U4  -  U3 

2  tan  9 

wi+w2  +  w3  +  w4  (u2  -  ui)  +  (y4  -  v3) 


(«3  +  «4)  ~  («1  +  «2)  (U3  +  V4)  -  (f7j  +U2) 

+^3 - — - V2 - 4 - 

.  w2  —  W\  w4  —  w3.l  +  cos2  9 
Hfs— - W— j— )-^2T 


(3) 
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V 


.  _  (tgl  +  t02)  -  (u>3  +  W4)  _  (u2  -  Ux)  -  (u4  -  Vg)  tan  ^ 


(«1  +  «2)  -  («3  +  Ui)  (Vl  +  V2)  ~  (V3  +  V4) 

+¥>3 - : - <P2 - 


,  /  Wj-Wi 
+  (^3 - - f"  ^2 


4 

w4  - 103 


)  tan  6 


Clearly,  (u,v,w)  do  not  form  a  true  velocity  vector  in  general.  For  the 
special  case  of  a  velocity  field  that  is  instantaneously  homogeneous  in  the 
horizontal  plane  over  the  distances  separating  the  beams,  (3)  does  give  the 
true  velocity  vector  on  a  ping-by-ping  basis  and  e  is  identically  zero. 

In  a  turbulent  flow,  however,  the  velocity  fluctuations  at  the  four  beams 
are  not  instantaneously  identical.  Thus,  the  instantaneous  horizontal  and 
verticals  components  of  the  velocity  cannot  be  measured  with  the  ADCP.  If 
the  scale  of  variations  of  channel  geometry  are  very  large  compared  to  the 
separation  of  the  beams,  then  the  statistical  properties  of  the  flow  can  be 
assumed  to  be  horizontally  homogeneous.  In  particular,  we  will  assume  that 
the  mean  components  (horizontal  and  vertical)  are  horizontally  homogeneous 
over  the  spatial  domain  of  the  beams  and  that  the  fluctuations  average  to 
zero.  Thus,  we  decompose  the  velocity  components  (u,-,Ui,  W{)  at  the  four 
beams  into  a  mean  field  (u,v,w),  which  is  the  same  for  all  the  four  beams, 
and  the  turbulent  fluctuations  (u',  v'^w';),  which  all  average  to  zero  but  may 
be  different  at  the  four  beams.  Then,  by  averaging  out  the  turbulent  fluctu- 
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ations,  (3)  becomes 


u  —  u 

V  =  V 

w  =  w  (4) 

e  =  e  =  0 

Hence  the  mean  velocity  vector  ( u,v,w )  can  be  obtained  by  sufficient  av¬ 
eraging  under  the  assumption  of  statistical  homogeneity  in  the  horizontal 
plane.  The  amount  of  averaging  required  to  gives  satisfactory  estimates  of 
the  mean  components  is  difficult  to  predict  a  priori.  The  amount  required 
will  certainly  depend  upon  the  distribution  of  eddy  scales  and  the  energy  of 
the  eddies.  For  example,  if  all  eddies  have  scales  large  compared  to  the  beam 
separations,  then  the  field  is  instantaneously  homogeneous  and  the  statistical 
uncertainty  of  the  estimates  is  determined  solely  by  the  probability  distri¬ 
bution  of  the  fluctuations  and  the  degrees  of  freedom.  For  such  eddies,  the 
vertical  velocity  must  be  small  compared  to  the  horizontal  components  owing 
to  its  bottom  boundary  condition.  The  most  intense  vertical  velocities  will 
come  from  eddies  with  scales  comparable  to  the  distance  from  the  bottom 
which  is  almost  the  same  as  the  distance  separating  b<=am  pairs.  The  aver¬ 
aging  will  have  to  span  many  such  eddies  which  is  realized  with  an  averaging 
time  t  L/U  where  L  is  the  beam  separation,  which  ranges  from  2.3  to 
31  m,  and  U  is  the  mean  speed.  Averaging  over  a  fixed  length  scale  is  not 
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very  practical.  We  chose  an  averaging  time  of  20  minutes.  The  least  amount 
of  spatial  averaging  occurs  for  weak  flows  and  amounts  to  only  120  m  in  a 
current  of  0.1  m  s-1.  For  typical  flows  of  1  m  s-1,  the  averaging  scale  is  1200 
m  and  the  ratio  of  averaging  scale  to  beam  separation  is  39  and  520  for  the 
top  and  bottom  of  the  profiles,  respectively. 

The  only  explicit  test  of  the  averaging  is  provided  by  the  redundant  es¬ 
timates  of  vertical  velocity  from  the  two  beam  pairs.  The  error  velocity  e  in 
(2)  and  (3)  is  half  the  difference  between  these  two  estimates  of  the  vertical 
velocity.  The  averaging  is  insufficient  or  the  assumption  of  horizontal  homo¬ 
geneity  is  invalid  if  the  error  velocity  is  not  small  compared  to  the  estimated 
vertical  velocity.  The  error  velocity  is  not  always  small  compared  to  the 
vertical  velocity  in  the  upper  half  of  the  profiling  range  during  weak  flows. 

For  a  rigidly  mounted  ADCP,  the  values  of  tilt  angles  <pi(i  =  1,2,3)  are 
constant  and  the  beam  velocities  can  be  averaged  first  to  estimate  the  mean 
velocity  components  from  (2),  that  is, 

—  &2  —  W  h  f  &2 

U  = - (£>3 - 

2  sin  8  2  cos  6 

—  64  —  63  63  +  64 

V  — - L.  - 

2  sin  6  2  cos  6 

—  61  +  62  +  63  +  64  62  —  61  64  —  63 

w  = - -7—=—a - +30  -  .  +<l2— —r  (5) 

4  cos  0  2  sin  V  2  sm  6 

|  _  (6j  +  62)  —  (63  +  64) 

4cos0 

Because  the  velocity  along  the  beams,  =  1  to  4),  can  be  measured  instan- 
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taneously,  the  uncertainty  of  their  means  can  be  evaluated  from  the  standard 
deviation  of  the  along-beam  fluctuations  h\.  The  statistics  of  the  along-beam 
means  can  then  be  used  to  estimate  the  confidence  interval  for  the  mean  flow 
vector  (u,v,w). 

Would  a  moored  ADCP  give  results  comparable  to  a  rigidly  mounted 
unit?  A  moored  instrument  would  have  fluctuating  tilt  angles  but  a  close 
examination  of  (3)  indicates  that  this  is  unlikely  to  have  significant  effects. 
All  products  involving  mean  tilt  angles  and  mean  velocities  would  vanish  just 
like  they  do  for  a  rigidly  mounted  instrument.  Products  involving  fluctuating 
tilts  and  mean  velocity,  and  vice  versa,  vanish  by  assumption.  Only  a  po¬ 
tential  correlation  between  fluctuating  tilt  angles  and  velocity  differentiates 
the  two  mounting  configurations.  For  example,  the  term  ^{u'2  —  u[)  from 
the  first  line  in  (3)  may  be  non-zero.  The  location  of  the  velocity  fluctua¬ 
tions  is  removed  from  the  position  of  the  instrument  and  only  fluctuations 
with  very  large  scales  are  likely  to  be  coherent  at  the  two  positions.  The 
coherence  would  be  equally  good  at  the  position  of  the  two  velocity  mea¬ 
surements  which  makes  the  correlation  small.  Similarly  if  u’2  and  u[  are  not 
coherent  then  neither  is  likely  to  be  coherent  with  tilt  fluctuations  which, 
again,  makes  the  correlation  of  tilt  and  velocity  small.  However,  the  tilting 
of  a  moored  instrument  is  unlikely  to  be  centered  around  its  transducers  and 
translational  movement  of  the  transducers  is  inevitable.  The  translational 
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motion,  which  is  not  monitored  by  the  instrument,  will  add  vectorially  to 
the  beam  velocities  and,  thus,  the  measured  beam  velocities  will  be  contam¬ 
inated.  Translational  motion  will  be  highly  correlated  with  fluctuations  of 
tilt.  The  vertical  velocity,  which  is  small  compare  to  the  horizontal  compo¬ 
nents,  is  unlikely  to  be  well  estimated  by  an  instrument  buffeted  by  turbulent 
fluctuations.  Without  reliable  estimates  of  vertical  velocity,  the  error  velocity 
is  also  unreliable.  Thus,  it  may  be  impossible  to  test  either  the  assumption 
of  horizontal  homogeneity  or  the  adequacy  of  the  averaging  for  a  moored 
ADCP  in  a  highly  turbulent  flow. 

4.  Data  analysis 

a.  Correction  of  tilting  angles 

Inaccuracy  in  the  pitch  and  roll  angles  measured  by  the  tilt  sensors  introduces 
errors  in  the  derivation  of  velocity  vector  from  beam  velocities.  This  affects 
the  calculation  of  the  vertical  velocity  most  because  it  is  weak  compared  to 
the  horizontal  components  and  is  easily  contaminated.  The  measured  values 
of  <p2  and  <£>3  are  0.00°  and  —3.00°  respectively  throughout  the  experiment. 
Using  these  two  values,  however,  the  calculated  vertical  velocity  at  mid-depth 
includes  a  tidally  varying  signal  (figure  3,  thinner  lines).  This  signal  is  not  the 
true  vertical  flow,  but  rather  the  contamination  by  the  horizontal  velocity. 
By  choosing  combinations  of  different  values  of  tp?  and  </?3,  we  found  that  this 
tidal  signal  can  be  minimized  by  changing  the  roll  angle  from  —3.0°  to  —2.0° 
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(thicker  lines,  figure  3).  Hence  the  values  of  <£>1,2,3  used  in  the  calculation  are 
(150.3°,  0.0°, -2.0°). 

b.  Measurement  uncertainties 

The  measured  beam  velocity  contains  an  uncertainty  due  to  the  Doppler 
noise.  To  determine  the  level  of  this  noise,  an  in  situ  test  was  conducted 
by  mounting  the  ADCP,  looking  downward,  on  a  boat  drifting  in  the  almost 
slack  water  of  Saanich  Inlet.  The  signal  associated  with  the  motion  of  the 
boat  was  removed  by  subtracting  a  second-order  polynomial  fit  from  each 
along-beam  velocity  profile.  The  remaining  signals  are  taken  as  Doppler 
noise.  Figure  4  shows  that  the  standard  deviation  of  the  beam  velocity  noise 
is  0.01-0.015  m  s-1  up  to  a  profiling  range  of  40  m,  for  4-ping  averaged 
ensembles  with  1  m  bin  size.  The  consecutive  ensembles  are  uncorrelated. 
Hence  by  averaging  N  ensembles  together  the  noise  standard  deviation  is 
reduced  by  a  factor  of  N 1/2.  For  20-minute  averages,  the  uncertainty  in 
beam  velocity  is  5-7  xl0~4  m  s-1. 

For  the  mean  velocity  estimates,  turbulent  fluctuations  introduce  another 
source  of  uncertainty.  Figure  5  shows  a  sample  'of  one-day  long,  4- ping  aver¬ 
aged  beam  velocity  at  mid-depth  from  beam  No.l,  the  beam  oriented  in  the 
downstream  direction  during  the  ebb.  The  flow  is  southward  during  the  ebb 
and  northward  during  the  flood.  A  low-pass  4-th  order  butterworth  filter, 
with  zero-lag  and  a  cut-off  period  of  20  minutes,  was  applied  to  separate 
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the  low-  from  high-frequency  variations.  The  intensity  of  fluctuations,  as 
seen  from  the  high-frequency  residue,  is  not  stationary  throughout  the  tidal 
cycle.  Figure  6  shows  the  standard  deviation  of  beam  velocity  fluctuations, 
calculated  over  20-minute  intervals,  for  beams  No.l  and  No.3  at  two  depths 
(beam  No.3  was  oriented  approximately  in  the  cross-channel  direction).  The 
time- variation  in  the  standard  deviation  is  obvious.  The  rms  velocity  fluctu¬ 
ation  is  about  0.05  m  s-1  at  peak  flow,  reaches  up  to  0.1  m  s-1  during  the 
flow  turning  at  day  24.05,  and  decreases  to  0.015  m  s-1  at  the  weak  flood 
between  day  24.6  and  24.8.  Hence  turbulence  dominates  the  high-frequency 
fluctuations  shown  in  figure  5,  and  is  only  comparable  to  the  Doppler  noise 
during  the  weak  flood. 

Most  statistical  analysis  requires  data  that  is  stationary,  or  homogeneous, 
in  the  time  domain  of  interests.  To  examine  the  stationarity  over  short  time 
intervals,  we  plot  two  sections  of  the  20-minute  long  time  series  of  the  velocity 
fluctuations,  one  taken  during  the  strong  ebb  starting  at  day  23.91  (figure  7), 
and  the  other  taken  during  the  flow  turning  starting  at  day  24.03  (figure  8). 
The  high-pass  velocity  fluctuations  plotted  are  along  beams  No.l  and  No.3, 
at  the  heights  of  3.6  m  and  15.6  m  respectively.  The  turbulent  fluctuations 
can  be  viewed  as  stationary  for  the  20-minute  interval  during  the  strong  ebb, 
but  the  fluctuations  during  flow  turning  are  not  stationary.  A  strong  event 
lasting  about  100  seconds,  with  changes  in  beam  velocity  exceeding  0.4  m  s_1, 
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happened  during  flow  turning.  In  fact,  the  large  values  of  the  beam  velocity 
standard  deviation  during  this  time  interval  (figure  6)  is  attributable  to  this 
single  event. 

Figure  9  shows  the  auto-correlation  functions  (acf)  against  the  time  lag 
r,  for  the  beam  velocity  fluctuations  shown  in  figure  7.  The  uncertainty 
level  of  the  acf  estimates  These  correlation  functions  are  not  like  that  for  a 
white  noise.  A  de-correlation  time  scale  is  defined  as  the  time  lag  when  an 
acf  is  reduced  to  the  upper  band  of  its  uncertainty  level,  which  is  estimated 
to  be  0.1.  For  this  20-minute- long  data,  the  de- correlation  time  is  about 
35  s  for  the  velocity  along  beam  No.  1  at  3.6  m  above  the  bottom  (figure 
9(c))  and  is  about  15  s  for  the  other  three  acf’s.  These  de- correlation  time 
scales  correspond  to  the  scales  of  the  dominant  turbulent  eddies.  Following 
Taylor’s  hypothesis  and  using  the  20-minute  mean  horizontal  velocity,  the  de- 
correlation  time  scales  found  above  imply  that  these  eddies  have  (horizontal) 
spatial  scales  ranging  between  15  and  30  m  during  this  strong  flow  interval. 
In  comparison,  the  acf’s  of  the  velocity  fluctuations  for  the  20-minute  interval 
starting  at  day  24.6  during  the  weak  flood  look  more  like  that  for  a  white  noise 
(figure  10).  This  may  simply  due  to  the  Doppler  noise  which  is  comparable 
to  turbulent  intensities  during  this  time  interval. 

Choosing  15  sec  as  a  typical  de-correlation  time  gives  80  independent  sam¬ 
ples  in  a  20-minute  interval.  The  standard  deviation  of  the  rapidly  sampled 
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beam  velocity  is  of  0(0.05)  m  s-1  typically  during  the  strong  flow  (figure  6). 
For  20-minute  averaging,  the  rms  velocity  fluctuation  is  reduced  to  0.0055 
m  s-1,  8-11  times  larger  that  due  to  Doppler  noise.  The  standard  deviation 
of  the  20-minute  mean  beam  velocity  is  0.0017  m  s-1  during  the  weak  flood 
between  day  24.6  and  24.8,  and,  still  dominates  that  of  the  Doppler  noise. 

Including  terms  representing  the  uncertainty  in  the  mean  beam  velocity 
in  equation  (5),  one  can  deduce  that  the  standard  deviations  of  the  horizontal 
and  vertical  velocity  components  are  \/2  times  larger  and  \/3  times  smaller 
respectively  than  that  of  the  beam  velocity.  Thus  the  uncertainty  due  to 
turbulence  in  the  20-minute  mean  velocities  is  about  0.008  m  s-1  for  (u,  v) 
and  0.003  m  s-1  for  w  typically  at  strong  flow,  and  the  level  of  uncertainty 
during  the  weak  flood  between  day  24.6  and  24.8  is  three  times  smaller. 

The  non-zero  acf’s  at  considerably  long  lags,  as  shown  in  figures  9  and 
10,  reflect  the  eddies  with  sizes  larger  than  those  of  the  dominant  ones.  For 
the  time  interval  during  the  flow  turning,  shown  in  figure  8,  the  distinct 
event  lasting  100  sec  dominates  the  eddies  with  shorter  time  scales.  Figure 
5  shows  that  a  smoothing  of  beam  velocity  at  20-minute  cutoff  period  not 
only  keeps  the  tidal  signal,  but  some  eddy  fluctuations  with  periods  longer 
than  20  minutes  as  well.  To  smooth  out  the  fluctuations  due  to  these  eddies 
requires  averaging  over  scales  larger  than  the  size  of  these  eddies.  However, 
choosing  a  longer  averaging  period  may  conflict  with  the  requirement  to  keep 
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the  variations  at  tidal  frequencies, 
c.  The  choice  of  averaging  period 

In  dealing  with  unstationary  tidal  flow,  the  choice  of  an  appropriate  averaging 
time  is  usually  hard  to  make.  This  is  because  the  “spectral  gap”  separating 
turbulence  from  the  tidal  fluctuations  is  “something  which  a  priori  may  not 
be  known  to  exist  or  physically  need  to  exist”,  as  noted  by  Gross  and  Nowell 
(1983).  A  compromise  is  always  needed  between  choosing  a  sufficiently  long 
averaging  time  to  reduce  the  uncertainties  by  turbulent  eddies  and  choosing  a 
shorter  one  to  keep  the  tidal  variations.  For  measurement  using  the  ADCP, 
there  is  another  constraint,  i.e.,  the  assumption  of  statistical  homogeneity 
may  also  require  sufficient  averaging.  The  20-minute  averaging  period  chosen 
in  this  analysis  reduces  the  rms  velocity  fluctuations  due  to  the  Doppler  noise 
and  turbulence  to  considerably  low  levels,  compared  to  the  mean  velocity 
magnitude.  The  change  of  the  tidal  flow  during  the  20-minute  intervals  is 
obtained  in  the  next  section  by  taking  a  harmonic  analysis  for  all  the  data. 
As  was  discussed  in  section  3,  20-minute  averaging  corresponds  to  averaging 
over  spatial  scales  much  larger  than  the  beam  separating  distance  during 
strong  flow.  Choosing  a  shorter  averaging  period  is  still  acceptable  at  peak 
flow,  but  would  be  too  short  when  the  flow  is  weak.  Whether  or  not  choosing 
a  shorter  averaging  time  changes  the  estimate  of  the  bottom  stress  from  the 
log-layer  fitting  is  addressed  in  the  companion  paper  by  Lueck  and  Lu  (1996). 
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5.  Results 


a.  Tidal  and  non-tidal  variations 

The  4.5-day  long  data  collected  using  mode  4  spanned  8  days  during  the  ex¬ 
periment.  Figure  11  shows  the  20-minute  mean  horizontal  flow  at  mid-depth 
and  the  harmonic  fit  composed  of  four  major  tidal  constituents  (M2,  S,2,  K\ 
and  Oi)  and  the  zero-frequency  residue.  The  flow  was  rotated  into  two  com¬ 
ponents  along  and  transverse  to  the  direction  pointing  345°  from  north,  which 
is  approximately  along  the  channel  axis.  The  harmonic  fit  explains  95%  of 
the  variance  in  the  observed  flow  at  mid-depth  and  91%  at  4.6  m  above  the 
bottom.  The  record  length  of  the  data  is  not  long  enough  to  obtain  accu¬ 
rate  harmonic  constants  for  all  the  constituents  included.  The  point  to  make 
here  is  that  the  major  signals  in  the  observed  channel  flow  are  tidal.  The 
change  of  flow  by  the  tidal  variations  over  the  averaging  time  intervals  cho¬ 
sen,  here  being  20-minutes,  reaches  up  to  0.15  m  s-1  at  mid-depth  during 
weak  flow  when  the  acceleration/deceleration  rates  are  largest.  Hence  the 
change  of  flow  due  to  tidal  variation  is  much  larger  than  the  uncertainties 
due  to  Doppler  noise  and  turbulence. 

The  difference  between  the  data  and  the  harmonic  fit,  is  also  shown  in  fig¬ 
ure  11.  Combining  this  difference  with  the  zero-frequency  residue,  obtained 
from  the  harmonic  fitting,  gives  the  non-tidal  variations  of  the  channel  flow. 
The  non-tidal  variations  can  be  attributed  to  the  large-scale  eddies  not  be- 
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Table  1:  The  cross-  and  along-channel  components  of  the  zero-frequency 
residual  flow  and  the  mean  magnitude  of  the  difference  between  the  observed 
flow  and  the  harmonic  fit,  in  [ms-1].  The  error  bars  of  the  mean  magnitude 
are  the  95%  confidence  intervals  obtained  using  a  bootstrap  method. 


z  [m] 

residue 

mean  magnitude  of  the  “difference” 

all 

ebb 

flood 

15.6 

0.090  ±  0.007 

0.078  ±  0.010 

4.6 

■BiiiMBijmi 

0.143  ±  0.009 

0.116  ±0.011 

ing  smoothed  out  by  the  20-minute  averaging,  and  the  secondary  circulation 
generated  mainly  by  the  curvature  forcing  of  the  headland.  Table  1  lists 
the  cross-  and  along-channel  components  of  the  zero-frequency  residue  and 
the  mean  magnitude  of  the  difference  between  the  observed  flow  and  the 
harmonic  fit  at  two  heights,  i.e.,  at  mid-depth  and  4.6  m  above  the  seabed 
respectively.  The  magnitude  of  the  “difference”  is  also  calculated  for  the 
ebb  and  the  flood  phases  separately.  The  error  bars  associated  with  the 
mean  magnitude  is  the  95%  confidence  interval  calculated  using  a  bootstrap 
method  (Efron  and  Tibshirani,  1993).  The  zero-frequency  residual  flow  is 
southward.  The  departure  of  the  flow  from  the  harmonic  fit  intensifies  near 
the  bottom,  and  is  larger  during  the  ebb  than  during  the  flood.  The  depar¬ 
tures  of  the  velocity  from  the  tidal  fit  are  attributed,  in  the  following  sections, 
to  asymmetry  of  the  secondary  circulation  with  respect  to  the  tide,  and  to 
large-scale  eddies  generated  by  flow  separation  from  the  headland  during  the 
ebb. 
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b.  Asymmetry  between  the  flood  and  the  ebb 

The  20-min  mean  velocity  profiles  of  (u(z),  v(z))  are  averaged  over  the  profil¬ 
ing  range,  and  referred  to  as  the  “depth-mean”  flow,  hereafter.  The  4.5-day 
long  data  of  the  depth-mean  flow  magnitude  and  direction  are  plotted  in  a 
polar  coordinate  diagram  (figure  12).  The  direction  of  the  ebb  is  approx¬ 
imately  160°  and  the  flood  is  approximately  350°,  i.e.,  the  two  directions 
differ  by  only  170°.  The  ebb  is  stronger  than  the  flood,  consistent  with  the 
net  southward  flow  obtained  from  the  harmonic  fitting.  The  current  direction 
also  changes  vertically  (figure  13).  The  direction  difference  between  bottom 
and  top  layers  is  about  15°  typically,  both  for  flood  and  ebb.  The  veering 
with  height  during  the  flood  is  in  the  same  direction  as  for  a  bottom  Ekman 
layer,  but  the  backing  during  the  ebb  is  in  the  opposite  sense.  The  vertical 
rotation  of  the  shear  (also  shown  in  figure  13)  is  about  90°,  much  larger  than 
for  the  flow.  Comparing  the  direction  of  the  shear  and  the  flow  at  3.6  m 
above  the  seabed,  one  can  see  that  both  are  quite  close  during  the  flood, 
whereas  the  shear  is  about  30°  counter-clockwise  from  the  flow  during  the 
ebb  (top  panel,  figure  13).  The  strong  asymmetry  in  the  direction  of  shear 
towards  bottom,  and  the  offset  of  shear  direction  from  the  current,  are  also 
illustrated  by  the  polar  coordinate  diagrams  of  the  4.5-day  long  data  (fig¬ 
ure  14).  The  significant  transverse  shear  during  the  ebb  corresponds  to  a 
transverse  bottom  stress  directed  offshore  from  the  headland  (figure  1). 
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c.  Velocity  and  shear  profiles 


The  direction  of  the  depth-mean  flow  is  not  constant.  Defining  the  stream- 
wise  direction,  s,  as  directed  downstream  of  the  depth-mean  flow  during  the 
flood,  and  the  transverse  direction  n  being  normal  and  to  the  right  of  s,  we  ob¬ 
tain  an  instantaneous  orthogonal  curvilinear  system.  Each  20-minute  mean 
velocity  profile  (u(z),  u(z))  is  rotated  to  components  (us(z),  un(z))  along  the  s 
and  n  axes.  The  streamwise  and  transverse  velocity  shear,  ( dus/dz ,  dun/dz ), 
are  calculated  from  the  (us(z),itn(z))  profiles  by  applying  a  first  difference 
for  the  top  and  bottom  bins  and  a  center  difference  for  the  middle  bins.  The 
uncertainty  due  to  turbulence  in  the  mean  shear  estimate  is  hence  0.008  s_1. 

Figure  15  shows  the  depth-time  sections  of  one-day  long  streamwise  and 
transverse  flow  («s,un)  and  shear  ( dus/ dz,dun/ dz ),  as  well  as  vertical  ve¬ 
locity  w.  Plots  of  time  series  and  profiles  of  these  quantities  are  provided 
(figures  16  to  19),  to  assist  the  interpretation  of  the  color  image.  The  stream- 
wise  flow  (top  panels  in  figures  15,  16)  shows  a  strong  tidal  signal,  with  the 
flow  in  the  lower  layer  diminished  by  the  bottom  drag.  The  non-tidal  vari¬ 
ations  can  be  strong  enough  to  distort  the  tidal  variations.  For  example, 
the  strong  ebb  between  day  24.3  and  24.6  has  two  peaks  and  with  weaker 
flow  in  between.  The  transverse  flow  is  stronger  in  the  near  surface  and  near 
bottom  layers,  and  weakest  at  mid-depth  (second  panels  in  figures  15,  16), 
Its  magnitude  reaches  0(0.15  m  s_1)  at  the  lowest  bin  during  the  ebb,  and  is 
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about  0(0.05  m  s-1)  in  both  the  upper  and  the  lower  layers  during  the  flood. 
The  transverse  flow  is  usually  directed  to  the  headland  in  the  lower  layer, 
and  in  the  opposite  direction  in  the  upper  layer,  with  exceptions  during  flow 
turning  and  the  weak  flood  between  days  24.6  and  24.8.  Also  during  the  ebb 
near  day  23.9  and  between  days  24.85  and  25,  the  transverse  flow  shows  a 
three-layer  structure.  Consecutive  evolution  (spaced  by  40  minutes)  of  veloc¬ 
ity  profiles  are  illustrated  in  figure  17.  The  maximum  speeds  of  streamwise 
flow  are  frequently  found  at  heights  below  the  upper  most  bin.  The  upper 
layer  apparently  leads  the  lower  layer  in  the  decelerating  stage,  whereas  it 
lags  during  the  accelerating  stage.  The  transverse  velocity  profiles  generally 
show  a  linear  variation  in  magnitude  with  height  in  the  lower  layer,  and  more 
complex  variations  in  the  upper  layer  (figure  17). 

Both  the  streamwise  and  the  transverse  shear  are  bottom-enhanced  (two 
lower  panels  in  figure  15;  figures  18  and  19).  The  streamwise  shear  reaches  a 
magnitude  of  0(0.05  s-1)  in  the  lower  layer,  and  contains  a  strong  tidal  signal. 
The  transverse  shear  reaches  about  half  the  magnitude  of  the  streamwise 
shear  near  the  bottom  during  the  ebb,  and  is  small  during  the  flood.  The 
streamwise  shear  frequently  reverses  sign  near  mid-depth,  in  accordance  with 
the  maximum  flow  speed  below  the  surface  (figure  17,  top  panel).  The  shear 
reversals  may  be  due  to  the  inertia  of  the  oscillatory  flow  (this  effect  on  near¬ 
bed  velocity  profile  was  discussed  by  Soulsby  and  Dyer,  1981).  However,  the 
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apparent  preference  of  shear  reversals  occuring  during  the  ebb  suggests  that 
they  could  result  from  the  entrainment  of  shallow  water  from  Saanichton  Bay, 
which  would  slow  down  the  water  above  mid-depth  in  the  channel.  Shear 
reversals  are  also  found  in  the  transverse  direction  during  the  ebb,  and  reflect 
the  occurrence  of  a  three-layer  structure  in  the  transverse  flow. 

d.  Vertical  velocity 

The  middle  panel  of  figure  15  shows  the  20-minute  mean  vertical  velocity. 
Several  significant  up-  and  down-welling  intervals  can  be  identified,  some 
happen  during  flow  turning  and  some  during  strong  currents.  These  inten¬ 
sified  vertical  flow  intervals  generally  occur  over  the  entire  profiling  range, 
although  the  maximum  velocity  is  at  mid-depth,  reaching  0(0.05  m  s-1)  at 
times.  Figure  20  compares  the  vertical  velocity  and  the  corresponding  “er¬ 
ror”  velocity  at  three  levels.  The  sign  of  the  “error”  velocity  is  kept  the 
same  as  that  of  the  vertical  velocity,  to  facilitate  the  comparison  of  these 
two  quantities.  The  magnitude  of  the  “error”  velocity  is  generally  smaller 
than  that  of  the  vertical  velocity  at  the  lowest  level.  This  indicates  that  for 
20-minute  averaging,  the  estimated  mean  velocity  is  homogeneous  over  the 
spread  of  the  beams  at  the  lower  levels  where  the  distance  separating  the 
beams  is  small.  The  magnitude  of  the  “error”  velocity  increases  with  the 
increasing  height.  At  the  first  peak  of  the  vertical  flow  after  day  24,  which 
corresponds  to  the  event  shown  in  figure  8,  the  “error”  velocity  is  large  and 
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comparable  with  the  vertical  velocity  at  15.6  and  27.6  m  above  the  bottom. 
Other  events  at  15.6  m  height  have  relatively  small  error  velocities  but  the 
error  velocity  is  frequently  large  at  27.6  m  above  the  bottom.  This  means 
that  20-minute  averaging  is  not  sufficient  to  give  satisfactory  estimates  of  th^ 
mean  flow  at  the  upper  levels  where  the  beam  spreading  is  large.  Note  an 
0(0.05  m  s-1)  vertical  flow  lasting  20  minutes  is  enough  to  bring  the  bottom 
water  to  the  surface  and  vice  versa. 

6.  Discussions 

a.  ADCP  deployment  techniques 

Various  techniques  to  deploy  the  ADCP  were  reported  among  the  increas¬ 
ing  applications  of  this  instrument  to  oceanic  measurements.  For  exam¬ 
ple,  ADCP  can  be  mounted  on  moving  ships,  looking  downward,  to  take 
transects  in  a  study  area  (e.g.,  Candela  et  al,  1992;  Geyer,  1993;  Gargett, 
1994),  mounted  in  a  buoy  (e.g.,  Schott  and  Johns,  1987;  Johns,  1988;  van 
Haren  et  al,  1994),  in  a  trawl-resistance  package  (Tang  and  Belliveau,  1994), 
or  mounted  rigidly  to  the  bottom  (Lohrmann  et  al,  1990;  and  this  work). 
The  data  recorded  can  be  the  beam  velocities  or  the  velocity  vectors  trans¬ 
formed  by  the  internal  processors.  Ensemble  averaging  can  be  performed 
before  recording.  For  self-contained  deployment,  the  battery  capacity  and 
the  memory  space  constrain  the  experiment  duration,  the  length  of  the  av¬ 
eraging  interval,  and  the  number  of  pings  that  can  be  averaged  into  an  en- 
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semble.  The  direct-reading  deployments  as  adopted  in  this  study  are  free 
from  these  constraints,  but  there  are  still  choices  to  be  made  in  setting  the 
profiling  parameters.  For  either  self-contained  or  direct- reading  deployments, 
the  profiling  parameters  set  should  ensure  that  the  recorded  data  can  derive 
meaningful  results. 

The  derivation  of  velocity  vectors  from  the  beam  velocities  is  based  on 
the  assumption  that  the  flow  field  is  homogeneous  over  the  horizontal  dis¬ 
tance  separating  the  beams.  In  a  turbulent  flow  field,  the  homogeneity  holds 
only  statistically.  Hence  even  though  the  beam  velocities  can  be  sampled 
at  rapid  rates,  only  the  time  averaged  velocity  vectors  can  be  derived.  For 
rigidly  mounted  deployment,  averaging  the  beam  velocity  to  derive  the  mean 
velocity  vector  gives  the  same  results  as  averaging  the  individually  trans¬ 
formed  “velocity  vector”.  For  non-rigid  mounting  deployment,  however,  the 
transformation  must  be  done  from  ping  to  ping  as  the  orientation  of  the 
beams  changes  with  time  hence  averaging  the  beam  velocity  is  meaningless. 
The  vertical  flow  is  usually  weak  and  is  easily  contaminated  by  small  er¬ 
rors  in  the  tilt  angles,  and  the  motions  of  a  platform  to  which  the  ADCP  is 
mounted.  For  a  rigid  mount  without  platform  motions,  a  systematic  bias  in 
the  measurement  of  tilt  angle  ^.an  be  corrected  through  minimizing  the  mean 
vertical  flow  signal  by  changing  the  values  of  tilt  angles.  The  vertical  velocity 
is  difficult  to  obtain  for  non-rigid  deployments  as  the  motions  of  the  unstable 
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platforms  such  as  ships  or  buoys  are  usually  not  easily  determined.  Thus,  it 
may  not  be  possible  to  test  the  homogeneity  of  the  mean  flow  estimates  by 
comparing  the  error  velocity  against  the  vertical  velocity.  The  measurement 
of  vertical  flow  using  a  modified  ADCT  with  a  true  vertical  beam  (Gargett, 
1994)  also  requires  a  sufficiently  stable  platform. 

b.  Channel  flow  dynamics 


The  measurement  revealed  a  complicated  picture  of  tidal  flow  in  a  natural 
channel.  The  curvature  of  the  side  boundaries  of  the  channel,  seen  in  figure 
1,  contributes  to  this  complexity.  Taking  account  of  flow  curvature,  the 
momentum  balance  in  the  transverse  flow  direction, n,  is 

'2  dy  8 


u: 


3u  du 

IT  + u-~d7  ~  R  -  =  ~9¥n  +  TC<w  ) 


(6) 


where  y  is  the  free-surface  elevation,  g  the  acceleration  due  to  gravity,  /  the 
Coriolis  parameter,  and  R  is  the  local  radius  of  streamline  curvature  (for 
depth- mean  flow).  The  effect  due  to  stratification  has  been  neglected.  The 
other  approximations  made  in  writing  (6)  can  mostly  be  justified  by  assuming 
0(un)  <  0(us )  (see  Kalkwijk  and  Booij,  1986).  The  second  term  in  (6), 
representing  the  advection  by  the  streamwise  flow,  was  not  measured.  The 
relative  importance  of  the  curvature  effect  to  Coriolis  force  can  be  quantified 
by  a  Rossby  number  R0  —  us/(fR).  For  /  ~  10-4  s-1  and  us  ~  1  m  s_1, 
R  —  2, 10  km  gives  R0  ~  5  and  1  respectively.  During  the  ebb,  un  generated 


26 


by  flow  curvature  is  in  the  direction  of  the  observed  secondary  circulation,  and 
it  is  opposite  to  the  contribution  by  Coriolis  force.  The  secondary  circulation 
generated  by  Coriolis  force  and  curvature  effect  are  along  the  same  direction 
during  the  flood,  and  consistent  with  the  observed  circulation.  The  stronger 
transverse  circulation  observed  during  the  ebb  suggests  that  the  streamlines 
have  sharper  curvature  during  this  phase  of  the  tide  than  during  the  flood. 

The  radius  of  streamline  curvature  is  not  known  because  the  measure¬ 


ments  were  taken  at  a  fixed  site.  However,  a  tentative  guess  can  be  made 
using  the  formulae  derived  by  Kalkwijk  and  Booij  (1986).  Assuming  that 
the  log-layer  in  the  streamwise  direction  extends  to  the  surface,  they  got 


:  =  sign(us)^fc(^a) 
Unb  2  K2Rfb(h’a) 


where  unc  denotes  the  contribution  from  Coriolis  forcing  and  uni  from  the 

1  /2 

effect  of  curvature  (bending);  us  is. the  depth-mean  of  us]  and  a.  =  Cjf  //c. 


The  two  profile  functions,  fc  and  fb ,  are  nearly  linear  functions  of  z  except 
very  close  to  bottom  (their  figure  2).  Using  our  estimates  of  Cd  x  10-3 
from  fitting  log-layers  to  the  streamwise  velocity  profiles  (Lueck  and  Lu, 


1996)  gives  a  ~  0.15,  the  magnitude  of  fc  and  fb  in  both  the  near-surface 
and  bottom  layers  are  0(0.5).  In  Cordova  Channel  h  ~  30  m  and  /  10-4 

s-1,  for  us  ~  1  m  s-1  we  have  unc  ~  0.01  m  s-1.  Hence  the  Coriolis  force 


does  not  dominate  during  either  the  flood  or  the  ebb.  Using  the  observed 
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secondary  circulation  during  the  flood  and  the  ebb,  unb  rj  0.05,0.15  m  s-1 
correspond  to  R  «  4, 1.5  km  respectively.  This  implies  that  the  streamline 
curvature  may  have  a  radius  as  small  as  1.5  km  during  the  ebb  and  more  than 
4  km  during  the  flood.  These  scales  are  not  inconsistent  with  the  geometric 
and  topographic  features  of  Cordova  Channel.  During  the  ebb,  the  flow  from 
the  Saanichton  Bay  experiences  a  sharper  curvature  when  turning  around 
the  headland.  The  channel  is  fairly  straight  to  the  south  of  the  experiment 
site,  hence  during  the  flood,  the  flow  aligns  with  the  axis  of  the  Channel  and 
can  over  shoot  the  headland  which  produces  a  smaller  curvature. 

Geyer  and  Signell  (1991),  also  Geyer  (1993),  observed  flow  separation  for 
tidal  flow  passing  a  headland.  The  generation  of  horizontal  eddies  by  tidal 
flow  was  also  observed  by  Wolanski  et  al  (1984)  in  the  lee  of  an  island.  The 
sharp  curvature  experienced  by  the  flow  during  the  ebb  in  Cordova  Channel 
may  also  favor  the  separation  of  flow,  as  suggested  by  visual  observations  in 
the  Channel  (Eric  D’Asaro,  personal  communication,  1995). 

7.  Summary 

In  this  study,  a  broadband  ADCP  was  rigidly  mounted  to  the  bottom  of 
Cordova  Channel  and  the  data  were  directly  read  by  a  computer  on  the 
shore.  The  rapid  sampling  of  beam  velocity  profiles  allows  for  an  explicit 
estimate  of  the  level  of  the  uncertainty  due  to  turbulence  in  the  calculated 
mean  velocity.  It  is  argued  that  for  flow  measurement  in  highly  turbulent 
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environment,  averaging  over  sufficient  large  scales  is  required  in  accordance 
with  the  assumption  of  the  flow  being  statistically  homogeneous  over  the 
distance  separating  the  beams  of  the  AD  CP.  The  20-minute  averaging  period 
chosen  in  this  analysis  corresponds  to  averaging  over  scales  large  compared 
to  the  beam  spreading,  and  also  reduces  the  uncertainties  in  the  mean  flow 
and  shear  estimates  caused  by  the  Doppler  noise  and  the  dominant  turbulent 
eddies.  The  stable  platform  provided  by  the  rigid  mounting  enables  removing 
the  bias  in  the  vertical  velocity  due  to  a  systematic  error  in  the  measurement 
of  tilt  angles. 

The  depth-time  variations  of  flow  and  shear  observed  in  Cordova  Channel 
are  complicated.  A  harmonic  analysis  including  the  major  tidal  constituents 
explains  more  than  90%  of  the  variance  in  the  20-minute  mean  velocities,  indi¬ 
cating  that  the  flow  in  the  Channel  is  mainly  tidal.  Separating  the  horizontal 
flow  into  streamwise  and  transverse  components,  substantial  transverse  flow 
and  velocity  shear  were  observed,  showing  a  strong  influence  of  flow  curvature 
provided  by  the  headland  at  one  side  of  the  channel.  The  curvature  effect 
dominates  the  Coriolis  forcing  in  driving  the  transverse  circulation.  The 
streamwise  shear  is  bottom-enhanced.  The  shear  frequently  reverses  sign  at 
mid-depth  during  the  ebb,  and  this  is  attributed  to  the  entrainment  of  the 
water  from  a  shallow  bay  to  the  north  of  the  channel.  Intensified  up-  and 
down-welling  events  with  magnitude  reaching  0.05  m  s-1  at  mid-depth,  are 
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observed.  These  large  values  of  vertical  velocity  are  validated  by  comparing 
them  with  the  corresponding  “error”  velocities. 
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Figure  captions 

Figure  1.  Area  map  and  topography  (depth  in  meters)  of  Cordova  Chan¬ 
nel  and  the  location  of  ADCP.  The  positions  of  two  deployments  of 
the  moored  microstructure  instrument  (TAMI1  and  TAMI2),  a  current 
meter  (CMI),  and  the  transmitter  (T)  and  receiver  (R)  of  a  acoustic 
scintillation  system  deployed  during  the  experiment  are  also  shown. 

Figure  2.  The  probing  definition  of  ADCP.  6(=  30°)  is  the  beam  angle  and 
2  and  ip3  are  heading,  pitch  and  roll  angles. 

Figure  3.  One-day  long  vertical  velocity  at  mid-depth,  calculated  for  (roll 
angle)  ip3  =  —2°  (thicker  lines)  and  <p3  =  —3°  (thinner  lines).  The  two 
upper  curves  are  calculated  from  20-minute  smoothed  beam  velocities 
and  the  lower  curves  are  3-hour  smoothed.  The  lower  panel  is  a  stick 
diagram  of  the  horizontal  flow  at  mid-depth. 

Figure  4.  Standard  deviation  of  the  4-ping  averaged,  1-m  bin  size  of  veloci¬ 
ties  along  the  four-beams,  and  its  variation  with  the  distance  from  the 
transducer.  ADCP  was  mounted  on  a  boat  drifting  in  the  almost  slack 
water  of  the  Saanich  Inlet  and  the  measurement  was  taken  using  the 
mode  4. 

Figure  5.  A  sample  of  one-day  of  velocity  data  at  mid-depth  from  the  beam 
oriented  in  the  downstream  direction  during  ebb.  The  sampling  time 
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interval  is  3.05  s  and  each  sample  is  the  average  of  4  pings.  The  20- 
minute  smoothed  velocity  and  the  residual  high-frequency  component 
are  plotted.  The  stick  digram  in  the  lower  panel  is  the  20-minute  mean 
horizontal  velocity  at  that  depth. 

Figure  6.  Time  variation  of  the  standard  deviation  along  beams  No.l 
(thicker  lines)  and  No.3  (thinner  lines),  at  mid-depth  (top  panel)  and 
3.6  m  above  the  sea  bed  (middle  panel).  The  bottom  panel  is  a  stick 
diagram  of  the  horizontal  flow  at  mid-depth. 

Figure  7.  Time  series  of  the  high-pass  velocity  fluctuations  along  along 
beams  No.l  (top  panel)  and  No.3  (bottom  panel),  at  mid-depth  (thin¬ 
ner  lines)  and  3.6  m  above  the  seabed  (thicker  lines).  This  is  a  20- 
minute  section  during  the  strong  ebb  starting  at  day  23.91. 

Figure  8.  Same  with  figure  7  but  during  the  turning  of  flow  starting  at  day 
24.03. 

Figure  9.  Auto-correlation  coefficients  of  the  four  high-pass  beam  velocity 
time  series  shown  in  figure  7. 

Figure  10.  Same  u.s  figure  9  but  for  the  20-minute  data  starting  at  day  24.6 
during  the  weak  flood. 

Figure  11.  The  cross-channel  (upper  panel)  and  along-channel  (lower  panel) 
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velocity  components  measured  at  mid-depth  (dots),  and  their  harmonic 
fits  including  M2,  S2,  Ki,  0\  and  a  zero-frequency  residue  (solid  curves). 
Differences  between  the  measured  flow  and  the  harmonic  fits  are  also 
shown  by  dots  with  an  offset  of  -1.5  m  s-1.  The  cross-  and  along-channel 
directions  are  pointing  to  75°  345°  from  the  north  respectively. 

Figure  12.  A  polar  coordinate  diagram  of  4.5  days  of  depth-averaged  20- 
minute  mean  velocity. 

Figure  13.  Upper  panel:  Direction  of  depth-mean  velocity  (thicker  solid 
line)  v.  flow  at  3.6  m  above  seabed  (thinner  solid  line).  Middle  panels: 
Typical  depth-variation  in  flow  direction  during  flood  (left)  and  ebb 
(right).  Bottom  panel:  The  change  in  shear  direction  with  depth. 

Figure  14.  Polar  coordinate  diagrams  of  4.5-day  long  flow  (left)  and  shear 
(fight)  at  Z  =  3.6  m  height. 

Figure  15.  Color-coded  depth-time  sections  of  (a)  streamwise  and  (b)  trans-  . 
verse  velocities;  (c)  20-minute  mean  vertical  flow;  (d)  streamwise  flow 
shear  with  the  log-layer  height  (dark  curve);  and  (e)  transverse  velocity 
shear. 

Figure  16.  Time  series  of  streamwise  (upper  panel)  and  transverse  flow 
(lower  panel)  at  the  heights  of  3.6  m  (thicker  solid  lines),  15.6  m  (thin¬ 
ner  solid  with  dots)  and  27.6  m  (thinner  solid). 
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Figure  17.  Profiles  of  streamwise  (upper  panel)  and  transverse  flow  (lower 
panel).  The  number  of  each  profile  indicates  the  index  in  the  sequence 
of  20-minute  segments  starting  from  day  23.91. 

Figure  18.  Same  with  figure  16  but  for  flow  shear. 

Figure  19.  Same  with  figure  17  but  for  flow  shear. 

Figure  20.  Upper  panel:  20-minute  mean  vertical  (thicker  lines)  and  “error” 
velocities  (thinner  lines)  at  the  heights  of  15.6m,  3.6  m  (offset  by  -0.05 
m  s-1)  and  27.6  m  (offset  by  0.05  m  s-1).  The  sign  of  the  ‘error” 
velocity  is  changed  the  same  with  the  vertical  velocity.  Lower  panel: 
Stick  digram  of  20-minute  mean  horizontal  velocity. 
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An  acoustic  Doppler  current  profiler  (ADCP)  was  applied  to  measure  the  three-dimensional 
flow  field  in  a  tidal  channel  along  the  coast  of  Vancouver  Island,  British  Columbia.  The  ADCP 
was  rigidly  mounted  to  the  bottom  in  the  center  of  the  channel,  and  rapidly  sampled  the  velocity 
profiles  along  its  inclined  beams.  The  data  were  directly  read  by  a  computer  on  shore.  This  method 
of  deployment  allows  for  an  explicit  estimate  of  the  uncertainty  of  mean  velocity  due  to  turbulence, 
and  obtains  the  weak  signal  of  vertical  velocity.  The  mean  flow  vector  is  derived  from  the 
measured  beam  velocities  under  the  assumption  that  the  flow  is  statistically  homogeneous  in  the 
horizontal  plane  over  the  distance  separating  the  inclined  beams  of  the  ADCP. 

The  depth-time  variations  of  velocity  and  shear  in  this  natural  channel  are  more  complex 
than  is  found  in  idealized  one-dimensional  channel  flow,  which  seldom  occurs  in  nature.  The 
channel  contains  strong  secondary  circulation,  intense  up-  and  down- welling  intervals,  frequent 
shear  reversals  and  a  substantial  amount  of  transverse  velocity  shear. 
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